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Multiple parton interactions in a single proton-proton collision are expected to play an important 
role for many observables at LHC. To a large part their phenomenological description relies on 
rather simple and physically intuitive assumptions. We report on an investigation of multiple hard 
interactions in QCD, which aims at identifying to which extent this simple description arises from 
theory and where it needs to be extended. An approximate connection with generalized parton 
distributions may help elucidate some aspects of multiple-interaction dynamics. 
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1. Introduction 

In proton-proton collisions at very high energies one can have events in which several par- 
tons in one proton scatter on partons in the other proton and produce particles of large transverse 
momentum or large mass. The effects of such multiparton interactions average out in sufficiently 
inclusive observables, which can be described by conventional factorization formulae that involve 
a single hard scattering. Multiple interactions do however change the structure of the final state and 
may be important for many analyses at LHC |jl|]. A brief review on the subject can be found in |Q]. 

Phenomenological estimates of multiparton dynamics rely on models that are physically in- 
tuitive but rather simplified, whereas a systematic description in QCD remains a challenge. The 
present contribution reports on some steps towards this goal. The discussion is limited to tree-level 
graphs, which already exhibit important features. Loop corrections, in particular soft gluon ex- 
change, add a further layer of complexity and are essential for determining if and in which form 
factorization holds. For details and further discussion we refer to ^. 

2. Cross section formula and multiparton distributions 

An example process where multiparton interactions contribute is the production of two elec- 
troweak gauge bosons (W, Z or y* with large virtuality) in kinematics where the transverse mo- 
menta of the bosons are small compared with their mass or virtuality. Since we are interested in the 
details of the final state, we keep the cross section differential in the transverse boson momenta. For 
the corresponding process with a single boson, a powerful theory in terms of transverse-momentum 
dependent parton densities has been developed |Q], which one may hope to generalize to the case 
of multiple hard scattering. 

Figure |I|a shows a graph with two hard-scattering subprocesses. The transverse momentum 
of each produced boson is the sum of the transverse momenta of a quark and an antiquark. With 
two partons emitted by each proton the transverse parton momenta are in general not equal in the 
scattering amplitude and its conjugate, as illustrated in the figure. Momentum conservation implies 
r + f = for the differences between the parton momenta in the amplitude and its conjugate. 
Performing a Fourier transform w.r.t. r and f one finds that their Fourier conjugate variables are 
equal, y = y. 




(b) 



Figure 1: Graphs for the production of two gauge bosons by double (a) or single (b) hard scattering. 
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Making the usual kinematic approximations in the hard-scattering subprocesses (i.e. neglect- 
ing small momentum components compared with large ones) one obtains 
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with 5 = 2 if the produced bosons are identical and 5 = 1 if they are not. a^^a, denotes the hard- 
scattering cross section for single-boson production, and the kinematic variables are s = {p + p)^, 
Xj = {qip)/{pp) andx,- = {qip)/{pp)- The definition of double-parton distributions closely resem- 
bles the one for the transverse-momentum dependent distribution of a single quark [Q]. For the 
two-quark distribution in the proton with momentum p we have 
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where we use light-cone coordinates m^ = (m° it u^)/^/2 for each four- vector u. One may regard k, 
as the "average" transverse momentum of each quark and y as their "average" transverse distance 
from each other, where the "average" refers to the scattering amplitude and its conjugate. Re- 
garding the transverse coordinates F(;c;,k,-,y) has the structure of a Wigner distribution [||], which 
depends on both momentum and position variables for each particle. (This does not contradict the 
uncertainty principle since Wigner distributions are not probability densities.) A similar application 
of this concept has been discussed for generalized parton distributions in [H]. 



Integrating (2.1) over the transverse boson momenta qi and q2 one obtains 
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The transverse-momentum integrated distribution /^([^^(.^oy) = S d^^\ ^^^2 Fai,a2{xijkhy) m^y 
be interpreted as the probability density for finding two quarks with momentum fractions xi and X2 
at a relative transverse distance y in the proton. The form (23) has long been known (see e.g. [^, |8]]) 
and underlies most phenomenological models for multiparton interactions. 



For each quark there are three relevant Dirac matrices in (2.2), 



r, = i7+, rA, = ir+r5, Pg^ = \iaJ+r5 with; = 1,2, (2.4) 



which respectively project on unpolarized, longitudinally polarized and transversely polarized 
quarks. Note that polarized two-parton distributions appear even in an unpolarized proton, since 
they describe spin correlations between the two partons. For small but comparable xi and X2 one 
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may well have sizeable spin correlations between the two quarks (which are close in rapidity for 
x\^X2), even if there is little correlation between the polarizations of a quark and the proton (which 
are far apart in rapidity). The relevance of such correlations in multiparton interactions was pointed 
out already in [||] but has to our knowledge never been included in phenomenological estimates. 

Note that if parton spin correlations are sizeable they can have a strong impact on observ- 
ables. For the production of two gauge bosons one can easily see that the product F^q^Aq pAqMi 
of longitudinal spin correlations enters the cross section with the same weight as the unpolarized 
term FqqFq^q. One also finds that product F8q,8q Fgq.sq of transverse spin correlations give rise to 
a cos(2<p) modulation in the angle (p between the decay planes of the two bosons and thus affects 
the distribution of final-state particles. 

The formulae given so far have ignored the color structure of the multiparton distributions. The 
quark lines with momenta k\ it ^r in Fig. |l|a can couple to a color singlet (as in single-parton distri- 
butions) but they can also couple to a color octet, provided that the lines with momenta k2±^r are 
coupled to a color-octet as well. Such color-octet distributions contribute to the multiple-scattering 
cross section, as was already pointed out in [||]. A more detailed discussion will be given in [||]. 



3. Power behavior 



It is easy to determine the power behavior of the cross section formula (2.1) for two hard 
scatters. The hard-scattering cross sections a behave like l/Q^, where 2^ ~ ^i ~ Q2 denotes the 
size of the large squared mass or virtuality of the gauge bosons. With transverse momenta qi ~ qi 
of generic hadronic size A we find 
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where we have used that the two-parton distributions scale like F ~ 1/A^ and that the typical 



transverse distance y between the partons is of order 1 /A. The same power behavior as in (3.1) 
is obtained for the case where both bosons are produced in a single hard scattering, as shown in 
Fig. |I|b. For the cross section differential in the transverse boson momenta, multiple hard interac- 
tions are therefore not power suppressed. 

The situation changes when one integrates over qi and q2. In the double-scattering mechanism 
both transverse momenta are restricted to be of size A, but for a single hard scattering one has 
|qi +q2| ~ A whereas the individual transverse momenta can be as large as Q. Because of this 
phase space effect one has 
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In the transverse-momentum integrated cross section, multiple hard scattering is therefore only 
a power correction. This is required for the validity of the usual factorization formulae, which 
describe only the single-scattering contribution. 
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4. Connection with generalized parton distributions 

A simple ansatz for modeling two-parton distributions is to write them as a product of single- 
parton distributions, thus neglecting correlations between the two partons. This provides a starting 
point for phenomenology, even though one may not expect such an approximation to be very pre- 
cise. A way to implement this ansatz for the distributions f (;c,-,k,-,y) is to insert a sum Y,x 1^) (^1 
over all physical states between the two bilinear operators in ( |2.2| ). If one assumes that the proton 
state is dominant in this sum, one has 

{p\ (^xj^.v,) {qxiqxi)\p) ~ J Y^Tl%\i {P\^^2qx2\p) {p\qx^qxi\p) , (4.1) 

where the subscripts x\ , X2 indicate the momentum fraction associated with each field and spinor 



indices have been omitted for brevity. After the Fourier transform in (2.2), momentum conser- 
vation imposes /?'+ = p^ but one still has p' ^ p. The two-parton distribution F(x;,k,-,y) is 
thus approximated by a product of two generalized parton distributions with zero skewness t, . 
For transverse-momentum integrated distributions one obtains a simple representation F(x,,y) = 
/ d^h /(x2 , b) /(xi , b + y) in terms of impact parameter dependent parton densities [^] . 

The same method can be applied to the color-octet distributions mentioned above. Rearranging 
the quark operators as 

\{qx2^''qx2) {qxi^^qxi) = - {qx,qx2) {qx2qx^)-\{qx2qx2) (qx^qxi) , (4.2) 

where X" are the Gell-Mann matrices, one can insert proton states between the operators in paren- 
theses on the r.h.s. For the first term this leads to generalized parton distributions with nonzero 
skewness t, , since the fields coupled to a color singlet carry different longitudinal momenta. 

Although it is only approximate, the connection with generalized parton distributions that can 



be measured in exclusive reactions — pointed out already in [10] — will hopefully help to better 
understand at least some features of multiple hard interactions. 
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